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Abstract 

(— I . Fixed-field Alternating Gradient (FFAG) accelerators offer the potential of high-quality, mod- 

I erate energy ion beams at low cost. Modeling of these structures is challenging with conventional 

beam tracking codes because of the large radial excursions of the beam and the significance of 
fringe field effects. Numerous tune resonances are crossed during the acceleration, which would 
lead to beam instability and loss in a storage ring. In a non-scaling FFAG, the hope is that these 
resonances can be crossed sufficiently rapidly to prevent beam loss. Simulations are required to sec 
if this is indeed the case. Here we simulate a non-scaling FFAG which accelerates protons from 31 
■ to 250 MeV. We assume only that the bending magnets have mid-plane symmetry, with specified 

vertical bending field in the mid-plane (y = 0). The magnetic field can be obtained everywhere 
^ I using a power series expansion, and we develop mathematical tools for calculating this expansion 

to arbitrary order when the longitudinal field profile is given by an Enge function. We compare 
the use of a conventional hard-edge fringe with a more accurate, soft-edge fringe field model. The 
tune 1/3 resonance is the strongest, and crossing it in the hard-edge fringe model results in a 21% 
loss of the beam. Using the soft-edge fringe model the beam loss is less than 6%. 

1 Introduction particles undergo from the lowest to the high- 

est energy, making the machines (and magnets) 
larger. 

The concept of a fixed-field, ahernating gradient jj^^^ interested instead in a non-scaling 

(FFAG) accelerator came to hght m the 1950's ppAG-in such a machine no attempt is made to 

[1, 2, 3], but has only recently been revived as a ^^^^ ^^^^^ constant during the acceleration. 

promising design for a low-cost, reliable acceler- jr.,-,. -j ui j • ^.i, 

^ ° ° ' in tact, the tunes vary considerably during the 

ator [4, 5]. Reasons for a lower cost are the fixed- acceleration, crossing many resonances [6, 7, 8]. 

field magnets (as opposed to those that must be rpi j. ■ i • i, • i j 

° ^ ine magnets m these devices can be simpler, and 

ramped up during each acceleration cycle), smah ^^^^.^^ excursion of the particles less. Simu- 
size, and only two different magnet types. ^^^-^^ ^^^^^ ^^^-^^^ ^^.^.^^1 inherently 

FFAG accelerators come in two varieties: seal- non-hnear. The success of the machine rests on 

ing and non-scaling. In a scaling FFAG the mag- the acceleration being rapid enough that the in- 

nets are designed so that the vertical and hori- stabilities do not have sufficient time to destabi- 

zontal tunes are constant over the acceleration lize the beam. 

phase. This guarantees stability of the particle Non-scaling FFAGs have been proposed for 

orbits, but at the expense of more complex mag- ^ ^^^^ accelerator [9, 10, 11] as weU as for 
nets which can be difficult to manufacture. An- 



other problem is the large radial excursion that 



hadron therapy [12, 13]. The world's first non- 
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Table 1: Design parameters for a sample FFAG. 



Element Type 1 (cm) Bq (T) G (T/m) 



DL 


drift 


40 






DS 


drift 


7.5 






BD 


cf bend 


22 


0.803952 


-12. < 


BF 


cf bend 


44 


0.555057 


8.0 



2 A sample FFAG 

We consider a simple FFAG for proton accel- 
eration consisting of 24 doublet sections. Each 
doublet consists of a pair of combined-function 
bend magnets ("cf bend") with parameters given 
in Table 1. We consider the magnetic field 
B(x,?/, z) = {bx,by,bz) where the vertical com- 
ponent by within the magnet has the form 

by{x,0,z) = {Bo + Gx)S{z), (1) 

in other words a vertical bending field Bq with a 
linear transverse gradient G. In the vertical co- 
ordinate y the magnet is assumed to have mid- 
plane symmetry. Mid-plane symmetry will be 
explained in detail in the next section, and we 
will generalize the transverse profile beyond lin- 
ear gradients. The function S{z) describes the 
longitudinal profile of the field. It is approxi- 
mately one inside the bending magnet, and de- 
creases smoothly to zero outside it. 

The transverse gradient G is large enough that 
the vertical bending field switches sign within 
the magnet. This is necessary to keep the in- 
crease in overall radius small over an order of 
magnitude increase in energy (Figure 1). 

These combined function magnets can also 
be considered as quadrupoles which have been 
shifted transversely. The field in the geometri- 
cal center of these magnets is not zero unlike a 
conventional quadrupole. In this sense they are 
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Figure 1: Top view of the FFAG ring, showing 
magnets and closed orbits from 31 to 250 MeV. 
The scale is in meters. 

a pair of focusing and defocusing quadrupoles. 
Figure 2 shows a single period of the lat- 
tice, with the small defocusing quadrupoles and 
larger focusing quadrupoles. These magnets are 
shown off center because they are conventional 
quadrupoles shifted off the beam-line. The dots 
in the centers of these magnets are the magnetic 
origins of these magnets, where the field is zero. 
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Figure 2: Top view of one period in the sam- 
ple FFAG ring, showing closed orbits from 31 to 
250 MeV. The scale is in meters. 

This design allows for a large energy accep- 
tance, from 31 MeV to 250 MeV protons, while 
the radius of the associated closed orbits ranges 



scahng FFAG, EMMA [14, 15] will allow scien- 
tists to study the unique beam dynamics of this 
kind of machine. 
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from 4 m to 4.2 m. Figure 1 shows closed orbits 
of 16 particles with energies ranging from 31 to 
250 MeV. 

Modeling this accelerator through the entire 
energy range is a challenge for any tracking code. 
The closed orbits range outward nearly 20 cm, 
so at the extremes of energy they are far from 
the center of the magnet. In Figures 1 and 2, 
fringe fields are modeled in the conventional way 
by giving each particle a position and momen- 
tum kick as it enters or exits a magnet. The 
position kicks contain terms proportional to the 
position of the particle cubed [19], which can 
become quite significant for a particle displaced 
0.2 m from the center of the magnet (Figure 4). 
In this paper we compare the effects of such a 
"hard-edge fringe" model with a more accurate 
field model where the rapid rise of the magnetic 
field in the fringe field region is directly inte- 
grated through. 




Figure 3: Horizontal entrance and exit angles 
(in degrees) versus momentum deviation for the 
defocusing quad (blue and red) , and the focusing 
quad (tan and green). 



Figure 3 shows the entrance and exit angles to 
the quads over the total range in energy. In the 
fringe field region, when these angles are large 
they provide significant vertical focusing. 
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Figure 4: Horizontal displacement in meters ver- 
sus momentum deviation, at the entrance and 
exit to the defocusing quad (blue and red), and 
the focusing quad (tan and green). Note that 
20 cm deviations occur. 




3 Magnets with 
symmetry 



mid-plane 



We consider a rectangular bend magnet in Carte- 
sian coordinates. The primary bending field is 
oriented along the y-axis, so that the bending 
takes place in the x-z plane. The primary bend- 
ing magnetic field B(rE, y, z) = {bx, by, bz) is Boy 
in the center of the bending magnet. This mag- 
net is assumed to have a specified vertical field 
in the mid-plane (y = 0), 



by{x,0,z) = T{x)S{z), 



(2) 



where T(x) is the transverse field profile and 
S{z) is a smooth function describing the longitu- 
dinal variation of the bending field through the 
magnet. We often will take the transverse field 
profile to be linear as in (1) but the formulas 
below give the general case. The longitudinal 
profile function S{z) is close to one within the 
magnet and decreases to zero outside it. S{z) 
can be determined by curve-fitting the field of a 
real magnet. A general function with this shape 
can be obtained using a ^-parameter Enge func- 
tion [20, 21], 



Siz;e) 



1 



1 + exp(2p(z/e)) 



(l-tanh(p(z/e))) (3) 
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where p{z) is a fc term polynomial of degree k—1, 

p{z) = oi + a2Z + a-sz'^ + . . . + akz''~^, (4) 

and e is a length scale over which the longitudinal 
field changes in the fringe region. Physically e is 
determined by the aperture of the magnet. Note 
that < S{z) < 1. We assume that k is even 
and that the highest order coefficient > 0, so 
that S{z) 1 as z ^ — oo and S{z) ^ as 

Z +00. 




-3-2-10123 



z 

Figure 5: Sample longitudinal profile functions 
S{z;e), withp(2;) = z (hyperbolic tangent, solid) 
as well as p{z) = (1 + 4a; — Sx"^ + x^)/A (dashed). 
The scaling factor e was set to 1. 

By adding two Enge functions, we can create 
a profile function of a magnet from z = a to b 

S{z) = -S{z - a;ea) + S{z - b;eb). (5) 

The locations a and b correspond to the physical 
entrance and exit of the magnet. Note that the 
hyperbolic tangent profile is a special case of the 
2-parameter Enge function with p{z) = z. Two 
Enge functions are shown in Figure 5, first the 
hyperbolic tangent and then a more general 4- 
parameter Enge function. 

The two profile functions S in (5) depend 
on the polynomial p{z) in (3), and even when 
€a = Cfe and the same polynomial p{z) is used, 
the field growth and decay will not be perfectly 
symmetric with respect to the center of the mag- 
net. This is due to the even powers in p{z). In 



order that the rise and fall of S{z) be symmetric 
we must change the sign of all the even terms 
between the entrance polynomial Pa{z) and the 
exit polynomial Pb{z)- The hyperbolic tangent 
polynomial p{z) = z contains no even terms and 
requires no modification. 

In order to provide sufficient room for the 
fringe field decay, we must increase the physical 
length of the "magnet" , at least a few units of e 
on each end. Use of the hyperbolic tangent func- 
tion {p{z) = z) preserves the integrated magnet 
strength, but the general Enge function does not 
and may require some adjustment to the overall 
strength Bq. 

We require the magnet to have mid-plane sym- 
metry, or that 

bx{x,y,z) = -bx{x,-y,z), (6) 
by{x,y,z) = by{x,-y,z), (7) 
bz{x,y,z) = -6^(x,-y,z). (8) 

If we consider T{x) and S{z) as known func- 
tions, then the symmetry constraints (6-8) plus 
Maxwell's equations determine the field every- 
where, and can be obtained using a power se- 
ries expansion. The magnetic potential ip can 
be written as the power series 

V' = yMx,z) + —fi{x,z) + —f2{x,z) + .. . (9) 

where B = VV' and V'^Tp = 0. Here the sym- 
metry dictates that ip must be an odd func- 
tion of y. The given transverse profile implies 
fo{x,z) = T{x)S{z) and in order that V'^ip = 0, 

ux,z) = (-i)"x; ('')t^''Hx)s^'-"''Hz). 

(10) 

In the special case with a linear gradient in the 
field, T{x) = Bq + Gx, this simplifies to 

/„(x,z) = (-l)"(i?o + Gx)S(2-)(z) (11) 

The power series (9) is valid for any x and 
z, but converges only when y is small enough. 
For the tanh profile function, p{z) = z, we can 
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obtain a bound on the radius of convergence. As 
shown in Appendix A, the field is guaranteed to 
converge for \y\ < R, where 

R=^e. (12) 

In the center of the fringe field region, R is equal 
to the radius of convergence for the power series 
(9). At any other value of x and z, it is a lower 
bound on the radius on convergence. 

If we allow the fringe fields to be arbitrarily 
short, e — 7- 0, the field diverges for arbitrarily 
small y. The reason for this divergence is that 
the rapid growth of the derivatives of the longi- 
tudinal profile function S^'^'^){z). It is not phys- 
ically realistic to allow the fringe field to decay 
over too short a length. The length of fringe 
field decay should instead be set according to 
the magnet aperture. 

Walstrom [22] has also computed the radius of 
convergence for magnetic field power series. He 
uses the on-axis generalized gradients, which are 
considered as complex- valued functions. The ra- 
dius of convergence at a particular location can 
then be determined by the nearest singularity in 
the complex plane. However, he does not calcu- 
late this radius of convergence for the case of a 
bend magnet with mid-plane symmetry. 

Venturini [23] has calculated higher order 
aberrations associated with hard-edge fringe 
field kicks in quadrupole magnets. He found that 
magnets with shorter fringe fields have larger 
high-order aberrations. In our view this can be 
explained by our calculation that the magnetic 
field power series diverges for all ?/ > as the 
fringe field length goes to zero. 

Instead of assuming a magnet with mid-plane 
symmetry, we could have used a conventional 
multipole field with the same longitudinal pro- 
file function S{z). The field is then a function of 
the radial coordinate r, the angle 6 and the lon- 
gitudinal variable z. We have repeated the same 
power series analysis, and find again that the 
magnetic potential power series is guaranteed to 
converge for |r| < R, where R is again given by 
(12). This multipole field may not converge for 



large values of the transverse variable x, whereas 
with the magnet with mid-plane symmetry the 
convergence of the power series (9) depends only 
upon y. Since transverse excursions of 0.2 m are 
present in our model FFAG accelerator, for soft- 
edge fringes we use field expansions with mid- 
plane symmetry rather than a quadrupole field. 

4 Modeling in PTC 

For modeling FFAGs, we use the Polymorphic 
Tracking Code (PTC), a beam tracking code de- 
veloped by Etienne Forest [24, 25]. For our pur- 
poses, the most important-indeed, essential — 
reasons for using PTC are its ability to place 
magnets in arbitrary locations (needed for the 
large transverse offsets of the quadrupoles), its 
correct handling of the large energy range in an 
FFAG, and its ability to integrate through arbi- 
trary magnetic fields. This last is done using an 
element of type arbitrary. 

Before any particle tracking, each arbitrary 
magnet is analyzed. Each component of the 
magnetic field in each transverse slice is fit by a 
bivariate polynomial ^ cijX^y^ with < < 
8. The fitting is done using a least-squares mini- 
mization. The magnetic field given by the poten- 
tial (9) is a polynomial in x and y, but not z. The 
fitting is necessary because the transverse slices 
are rotated from this orientation. To track a par- 
ticle through this element, PTC then integrates 
using these pre-calculated polynomial fields by 
means of a Lorentz kick. Here this integration 
scheme is not symplectic, but it is sufficient for 
short term tracking (< 1000 turns). 

In element of type arbitrary the fringe field re- 
gion is taken into account directly through the 
rise in the magnitude of the magnetic field, and 
not using a delta function kick. This is more 
accurate, but presents several difficulties. First, 
the integration step must be small in the fringe 
field, where the field is increasing rapidly. Sec- 
ond, the magnet must be longer to accommodate 
the fringe fields. The combined function bends 
are separated by a short drift of only 7.5 cm, so 
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that their fringe fields overlap. We have mod- 
eled this section by making the magnets 15 cm 
longer on each end, as can be seen in Figure 9. 
In order to patch the various overlapping com- 
ponents together, PTC automatically drifts the 
particles backwards to the beginning of the next 
component. These backward tracks can be seen 
in Figure 9. This method of integrating between 
overlapping components is justified by the fact 
that the field is small in the region of overlap. 

In all cases we use a field profile S{z) with a 
hyperbolic tangent profile, p{z) = z. The decay 
length e was set to 5 cm, which is similar to the 
aperture. We use an integration step of 1 cm 
within elements of type arbitrary. 

An alternate approach would be to consider 
the two bend magnets as one long magnet. A 
3D magnet code (such as TOSCA) can be used 
to generate the magnetic field of both magnets 
together as a single beamline element. We did 
not do this because it would entail additional 
complexity, for one thing the two magnets are 
rotated with respect to one another. Our model 
represents a compromise between a standard 
tracking code using hard-edge fringe kicks, and 
a model which tracks through many magnets us- 
ing the full 3D field. We retain the convenience 
of modeling the beamline simple sequence 
of elements (DL, BD, DS, BF) while reproducing 
more realistic fringe fields. 



5 Results 

We ran PTC both with the special field- 
mapped magnets, and also with conventional 
quadrupoles with a hard-edge fringe field en- 
trance and exit kick. These magnet types differ 
primarily in their treatment of the fringe field. 
Figure 6 shows the horizontal and vertical tunes 
for both magnet types. The tunes are similar 
but differ in key areas. In particular, the verti- 
cal and horizontal tunes are over 8 after injection 
for the conventional quadrupoles, but both are 
slightly under 8 for the field-mapped magnets. A 
machine tune of 8 corresponds to a single-period 
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Figure 6: Machine tune versus kinetic energy. 
The solid lines are for hard-edged magnets, the 
dashed lines for soft-edge (field mapped) mag- 
nets. 

tune of 1/3, because we have a 24 cell lattice. 
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Figure 7: Machine tunes versus kinetic energy. 
Solid lines unchanged, dashed lines are for soft- 
edge quads with increased strength G at 8.4 and 
-13.2 T/m. 

If we increase the strength of the quadrupoles, 
or equivalently increase the gradient of the field 
G, the tune diagrams become more similar to 
those for hard-edge magnets. Figure 7 shows the 
result when G is increased to 8.4 and -13.2 T/m, 
and Figure 8 when G is increased to 8.8 and -13.6 
T/m. For low energies the tunes are very close 
to the hard-edge case, and in particular the tune 
8 resonance must again be crossed. There are 
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still differences in the tunes at higher energies. 
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Figure 8: Machine tunes versus kinetic energy. 
Solid lines unchanged, dashed lines are for soft- 
edge quads with increased strength G at 8.8 and 
-13.6 T/m. 

Figure 9 shows closed orbits for both conven- 
tional hard-edge fringes and using the new soft- 
edge magnets. We can see that the horizontal 
focusing is decreased by the soft-edge fringes, so 
the transverse excursion of the orbits is greater. 

Figure 10 shows closed orbits for the magnets 
with the largest gradient (8.8 and -13.6 T/m), 
compared with the closed orbits of the original 
hard-edge fringes (same as in Figure 9). We can 
see that the horizontal focusing now appears sim- 
ilar, although the soft-edge orbits are displaced 
inwards, particularly at higher energies. 

We now take 1000 protons and accelerate them 
from 31 to 250 MeV using PTC. The start- 
ing emittance was 167r mm-mrad in each plane. 
The particles were uniformly distributed inside a 
phase space ellipse. In order to reduce the insta- 
bilities due to the many resonances that must be 
crossed, the acceleration needs to take place as 
rapidly as possible. For this reason we place an 
RF cavity in each of the long drifts. The 24 RF 
cavities add 350 KeV per turn, so acceleration 
takes 626 turns. The current in this machine is 
small enough that we do not take into account 
space charge. 

Using conventional hard-edge fringe field mag- 



nets, 21% of the beam is lost in the first 20 turns. 
This occurs at the crossing of the 1 /3 tune (Fig- 
ure 11). No other resonance gives further loss of 
the beam. Using the special field-mapped mag- 
nets, none of the beam is lost. This is not sur- 
prising because the tunes no longer cross the 1/3 
resonance. 
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Figure 11: Fraction of the beam lost by energy 
using hard-edge magnets. 

If we increase the quadrupole gradient of the 
magnets, the tunes again cross the 1/3 reso- 
nance, but a 21% beam loss is not seen using 
the soft-edge fringes. With the magnets of Fig- 
ure 7 only 3.4% of the particles are lost and with 
Figure 8 only 5.6% are lost. 

The hard-edge fringe field kick derived in [19] 
and used for quadrupoles in PTC assumes that 
the fringe field has quadrupole symmetry. A sub- 
tle point is that since our field is assumed to have 
mid-plane symmetry, the quadrupole fringe field 
kick may not be accurate for fringe fields with 
mid-plane symmetry. It may be that some of 
the particle loss differences between hard-edge 
and soft-edge magnets is due to this problem. 
Ideally, one would like to repeat the calculation 
in [19] for a magnet with mid-plane symmetry, 
and compare results using this fringe field kick 
to the soft-edge fringe model, but this is not a 
simple calculation and we have not done this. 
The fact that the closed orbits for hard and soft- 
edge magnets in Figure 9 are so close to one an- 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Figure 9: Top view of a single period showing closed orbits for hard-edge magnets (black, same as 
Figure 2) and soft-edge magnets (red) with G at 8.0 and -12.8 T/m. The scale is in meters. 



other suggests to us that our hard-edge fringe 
field kicks are reasonably accurate. 

Our results indicate that modeling fringe fields 
correctly is critical to predict the performance of 
non-scaling FFAGs. Although particle tracking 
with the new field-mapped, soft-edge magnets 
is about 10 times slower than using conventional 
hard-edge fringe model, the additional computa- 
tional time far outweighs any later re-engineering 
cost after a machine is built. In addition, the fact 
that we are only tracking particles through 626 
turns, rather than millions of turns as in a stor- 
age ring, means that the tracking time for 1000 
protons goes from seconds to minutes, in other 
words is still easily handled. 

6 Summary 

We have used the tracking code PTC to track 
protons through the acceleration phase of a non- 
scaling FFAG from 31 to 250 MeV. We have com- 
pared the use of a conventional, hard-edge fringe 
model with a more accurate, soft-edge fringe 
model. With an emittance of IGvr mm-mrad, the 
hard-edge fringe model predicts a 21% particle 
loss as the beam moves through the 1/3 tune 
resonance. In contrast the more accurate fringe 
field model predicts only a 6% particle loss. 



A future paper will study the effect of magnet 
misalignments for this machine [26]. 
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A Enge function derivatives 
and radius of convergence 

Differentiation of the Enge function (3) is 
straightforward but results in formulas of rapidly 
increasing complexity. The results are easier to 
comprehend if we interpret the Enge function as 
a composition of two functions, S{z; e) = f og{z) 
where 

f{z) = ^(1-tanhz) and (13) 
g{z) = p{z/e). (14) 

Derivatives of S can be found using Faa di 
Bruno 's formula [16] for differentiation of a com- 
position of two functions where the second sum 
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Figure 10: Top view of a single period showing closed orbits for hard-edge magnets (black, same 
as Figure 2) and soft-edge magnets (red) with G at 8.8 and -13.6 T/m. The scale is in meters. 



d" 



^/ o g{z) = ^ /C^) o g{z) J2 Hn; ai, a2, . . . , a„)(g'(z))"^ {g"{z)r ■ ■ ■ (5^"^ (^))°" (15) 



k=l 



part 



d" 



^/o5(^) = E/^'^°^(-)E 



n 



k=l 



ai lao! • • • a„! 

part ^ ^ i=l 



(16) 



in (15) or (16) is over the partitions a set of 
n = ai + 2o2 -|- . . . + nan different objects into 
subsets containing k objects for A; = 1, 2, . . . , n, 
and k = ai + a2 + ■ ■ ■ + an.. The coefficient 
$(n; ai, a2, . . . , a„) is the number of such par- 
titions. These are related to multinomial coeffi- 
cients and are defined by 



$(n;ai,a2, . . . ,0^) = JJ 



\ {il)''^a,\ 



(17) 



For (15) and (16) we need derivatives of 
tanh(^;) to arbitrary order. We find that 



d" tanh z 



Q'„(tanhz) 



where qn is a polynomial of degree n + 1 . These 
polynomials are defined by go = and the recur- 
rence relation 



In+liz) = {l-Z ) — 



(19) 



In [18], it is shown that these polynomials can 
be written in terms of Stirling Numbers of the 
second kind, 52 (n, A:). 



n 

qn{z) = {-2r{z + l)Y,S2{n,k)k\ 

k=0 



1 



(20) 



Table 2 gives the polynomials g„ up to n = 
12, and Figure 12 shows them graphically up to 
n = 10. Table 3 gives the coefficients needed in 
Faa di Bruno's formula (15) to calculate the 5th 
derivative of / o g{z). As an example, using the 
results of Tables 2 and 3, we calculate the 5th 
derivative of the Enge function S{z; e) as 



(18) d^S 



dz^ 



1 

2? 



p(5)gi + (5p(V + 10p(V)g2 

+ (10p(3VV + 15p"p"p')q3 
+I0p"{p'fqi + {p'fq^] 



(21) 



where derivatives of the polynomials p are eval- 
uated at z/e, and the polynomials q^ are evalu- 
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Table 2: The first 12 polynomials qn{z). 


n 


qn{z) 





Z 


1 




2 


-2^(1-^2) 


3 


-2(1 -z2)(l -3^2) 


4 


8z(l-z2)(2-3z2) 


5 


8(1 - z2)(2 - 15z2 + 15^4) 


6 


-16z(l - z2)(i7 _ 60z2 + 45z^) 


7 


-16(1 - z2)(i7 _ 231^2 + 525z^ - 315z^) 


8 


1282(1 - z2)(62 - 378^2 + 630z^ - 315^6) 


9 


128(1 - z2)(62 - 1320^2 + 5040z^ - OGlSz*' + 2835z^) 


10 


-256z(l - z2)(i382 - 12720^2 + 34965z^ - 37800^*^ + 14175z^) 


11 


-256(1 - z2)(i382 - 42306^2 + 238425^^ - 5093552*^ + 467775^^ - 1559252^°) 


12 


1024z(l - z2)(2l844 - 280731^2 + lUmOz'^ - 1954260^^ + 1559250z*^ - 467775z^°) 



0.5 - 



-0.5 





■\ / 






. . \ / 

y-., / 
^ ■■ \ 


\ J''^ ^ 
\ \ ■•■■■\-- 

/ V""/-!./..--' ' 
/ \/ 1 




vy \ \ 

\ '\ 

■■■■■■■■^ 


i / 



-1 



-0.5 



0.5 



Figure 12: The polynomials go(-z) to g'io(-z) nor- 
malized by i|-„/2]- 

ated at tanh(g((z)) = tanh(p(z/e)). In the spe- 
cial case of the tanh function, where -piz) = z, 
all the terms in (21) are zero except for the last 
one and we have 



d^tanh(z/e) g5(tanh(z/e)) 
dz^ 

which is the same as (18) for n = 5. 



(22) 



The absolute value of the lowest order co- 
efficient in qn{z) (the constant or linear term, 
whichever is nonzero) generates the sequence 1, 
1, 2, 2, 16, 16, 272, . . .. The sequence 1, 2, 16, 
272, 7936, 353792, ... is sequence A182 in OEIS 



Table 3: The seven terms in 5(5). 
k partition of 5 multin. coeff term 



1 

2 
2 
3 
3 

4 
5 



5 

4+1 

3+2 

3+1+1 

2+2+1 

2+1+1+1 

1+1+1+1+1 



$(5;0,0,0,0,l) 
$(5;1,0,0,1,0) 
$(5;0,1,1,0,0) 
$(5;2,0,1,0,0) 
$(5;1,2,0,0,0) 
$(5;3,1,0,0,0) 
^>(5;5,0,0,0,0) 



5p(4)p'g2 
10p(3)p"g2 



15(p")Vt 
10p"(p'^3, 
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[17], known as the "tangent numbers" tn- These 
are Taylor series coefficients of the tanh(2:) func- 
tion expanded about z = 0, and may also be 
written in terms of Bernoulli numbers [17]. In 
Figure 12 the polynomials qn{z) are divided by 

and it is clear that 



the normalizing factor t \n/2'] ) 

\qn{z)\ < for l^l < 1, 



(23) 



at least for n < 12. For odd n, (23) is equivalent 
to the statement that the maximum magnitude 
of the n'th derivative of tanh(z) occurs at z = 0. 
This can be proven by writing the Taylor expan- 
sion of the n'th derivative of tanh(z) about z = 
and applying the Alternating Series Estimation 
Theorem. 

Since the lowest order coefficient of quiz) is 
t|-„/2]) ill the center of the fringe field we can 
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say: 



References 



(0) = (-l)("-l)/2 / Hn+l)/2 n odd, 

" 10 n even, 



(24) 



and 



V«/2 





tn/2 



n odd. 



n even. 



(25) 



The asymptotic form for the tangent numbers 
is given by [17] 



2n 

tn = 2{^] (2n-l)!. (26) 



This formula is quite accurate even for small 
n — for example it predicts t4 = 271.958, when 
t4 = 272. We can use this asymptotic form to ob- 
tain a bound on the radius of convergence for the 
magnetic field power series (9) for the tanh lon- 
gitudinal profile and linear transverse field gra- 
dient (1). If the power series terms in (9) are c^, 
i.e. ip = "^CnU^, then = for even n and 
Equation (26) gives the bound 



< 



\Bo + Gx\ 
n{n — 1) 



vre 



forn > 1. (27) 



The radius of convergence of the power se- 
ries with coefficients c„ is given by \/R = 
limsup„_j.o^ Y^|cn|, which leads directly to the 
bound 



R 



limsup y\cn 

n— ^oo 



> 



vre 



(28) 



Thus the power series (9) converges for |y| < 
7re/2 (for any value of x and z). 
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